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APPLICATION OF THE AERODYNAMIC ENERGY CONCEPT 
TO FLUTTER SUPPRESSION AND GUST ALLEVIATION 
BY USE OF ACTIVE CONTROLS 

E. Nissim,* A. Caspi,^ and I. Lottati^ 

Langley Research Center 

SUMMARY 

The effects of active controls on flutter suppression and gust alleviation of two 
different types of subsonic aircraft (the Arava twin turboprop STOL transport and the 
Westwind twin-jet business transport) are investigated. The active control surfaces are 
introduced in pairs which include, in any chosen wing strip, a 20-percent chord leading- 
edge (L.E.) control and a 20-percent chord trailing-edge (T.E.) control. Each control 
surface is driven by a combined linear-rotational sensor system located on the activated 
strip. The control law is based on the concept of aerodynamic energy and utilizes pre- 
viously optimized control law parameters based on two-dimensional aerodynamic theory. 
The best locations of the activated system along the span of the wing are determined for 
bending-moment alleviation, reduction in fuselage accelerations, and flutter suppression. 
The effectiveness of the activated system over a wide range of maximum control deflec- 
tions is also determined. Two control laws are investigated. The first or basic control 
law utilizes both rigid-body and elastic contributions of the motion. The second or 
extended control law employs primarily the elastic contribution of the wing and leads to 
large increases in the activated control effectiveness as compared with the basic control 
law. 

The results indicate that flutter speed can be significantly increased (over 70 per- 
cent increase) and that the bending moment due to gust loading can be almost totally elim- 
inated by a control system of about 10 to 20 percent span with reasonable control-surface 
rotations. Large reductions in fuselage accelerations can also be obtained by a small 
number of activated L.E. -T.E. control systems. 


*NRC-NASA Resident Research Associate, on leave from Technion - Israel Institute 
of Technology, Haifa, Israel. 

^Israel Aircraft Industry, Lod, Israel, formerly with Technion — Israel Institute of 
Technology. Haifa, Israel. 

^Technion — Israel Institute of Technology, Haifa, Israel. 



INTRODUCTION 


The technological advances made in recent years in the field of control systems 
have stimulated considerable interest in evaluating the advantages of incorporating active 
control systems in aircraft. (See ref. 1.) The potentials of active controls for gust alle- 
viation and mode stabilization have been analyzed for some specific aircraft such as the 
XB-70 (refs. 2 and 3) and the B-52 (refs. 4 and 5). Within the last few years, control 
systems that increase the damping of the lower frequency structural modes have evolved 
from analytical feasibility studies to production hardware. Such a system, which controls 
the response of the rigid-body mode and one elastic mode (first aft body bending) to gust 
inputs, has been successfully installed on the B-52G and B-52H fleets. A reduction in 
gust loads and a considerable extension of the fatigue life of the aircraft (ref. 6) have 
resulted. Some recently developed hardware indicates that flutter -suppression systems 
(controlling higher frequency imstable modes) are now technologically feasible. Analyti- 
cal studies have shown that in many instances, weight savings by as much as 4 percent of 
the total structural weight of large aircraft like the supersonic transport or the Rockwell 
International B-1 can be achieved by suppressing flutter by active controls rather than by 
passive methods. This is a considerable weight saving when it is considered that the pay- 
load may be as small as 20 percent of the structural weight. (See refs. 7 and 8.) 

A number of investigators have recently embarked on experimental flutter- 
suppression studies aimed at verifying analytical results. These studies provided exper- 
imental verification regarding the possibilities of suppressing flutter by active controls 
(refs. 9 to 15) and led to the pioneering flight-test demonstration using the B-52 airplane. 
(See ref. 16.) 

One of the major difficulties which characterizes the introduction of active control 
systems into elastic structures lies in the need to determine a large number of param- 
eters associated with the control system. An additional difficulty is introduced by the 
fact that an elastic structure, like an aircraft, cannot be considered as a fixed system 
since the properties of the system vary with the flight configuration, time of flight, etc. 
Hence, a proper optimization process must sweep over a very large number of param- 
eters, including a large number of flight configurations. The determination of a satisfac- 
tory control law which copes with the variety of flight configurations can be seen to be a 
difficult task which requires considerable ingenuity from the control system designer. 

In an attempt to reach an optimization procedure which does not exhibit sensitivity 
to flight configurations and to reduce to a minimum the number of arbitrary assumed 
parameters, the aerodynamic energy concept was developed. (See refs. 17 to 19.) The 
following general results are obtained by using the aerodynamic energy analysis: 
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(a) Optimum systems should include both leading-edge (L.E.) and trailing-edge (T.E.) 
control surfaces in each controlled strip. 

(b) Optimum systems should be activated by both linear and rotational sensors 
located on the activated strip. 

(c) Very general control laws can be assumed (based on the nature of the aerody- 
namic forces), and their coefficients determined for optimum performance irrespective 
of the type of aircraft, mass and stiffness distribution, center of gravity (c.g.) or elastic - 
axis locations, and reduced frequency or Mach number (within the subsonic range) . 

The remaining parameters such as the spanwise location of the active strips and 
some free gain parameters can be determined for specific aircraft together with the mag- 
nitude of the improvement obtained. 

Flutter -suppression and gust-alleviation problems are closely interrelated and can- 
not be treated separately. This is true since flutter-suppression considerations impose 
no limitations on the values of the control law parameters whereas gust response consid- 
erations yield an upper bound to the control law values. This upper bound insures rea- 
sonable control- surface rotations over the flight envelope of the aircraft. Furthermore, 
the effects of the flutter-suppression system on the gust response characteristics of the 
aircraft need to be established. For this reason, the effects of the free parameters of 
this problem are investigated by bearing in mind both flutter -suppression and gust- 
alleviation problems. The purpose of the present work is to apply the initially optimized 
results, obtained through the use of the aerodynamic energy method, to flutter suppression 
and gust alleviation of specific aircraft. More specifically, the effects of the control 
parameters on the following aircraft characteristics are investigated: 

(1) The maximum absolute value of the bending moment acting on the wing due to 

gust input, 

(2) The acceleration at the center of gravity of the aircraft and along the fuselage 

due to gust input, and 

(3) The flutter speed. 

Some indication of the effects of the active control system on the short -period sta- 
bility is included in this paper. However, no consideration is given to problems associ- 
ated with controllability. 


SYMBOLS 


b semichord length 

[C],[G],[R] control law matrices defined by equation (1) 


3 



aerodynamic forces 
Fg control forces 

Fj inertia forces 

^RE ~ 

h displacement defined in sketch (a) 

h = h - hr 

bending moment 
t time 

V velocity 

Vg design cruising speed of aircraft 

Vg upgust velocity 

O', ft 6 rotations defined in sketch (a) 

a = Q! - Q!r 

6c value of 6 when either bending moment or center-of-gravity acceleration of 

unalleviated aircraft reaches maximum value 

w circular frequency 

c^rj^r fj^r 2 reference frequencies 

Subscripts: 

eg center of gravity 

f flutter 
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max parameter relates to maximum value 

r parameter relates to reference section 

Dots over S 3 mibols denote derivatives with respect to time. 

CONTROL SYSTEM 

The basic control system consists of a pair of leading- and trailing-edge control 
surfaces alined in the streamwise direction and having a width in the order of a tenth of a 
wing semispan. The control surfaces together with the wing surface between them are 
referred to as an activated strip. In the simplest form of the control concept, the deflec- 
tions of the control surfaces are commanded to be proportional to the linear and rotational 
motion of the strip. A control system consists of one or more activated strips located on 
the wing and/or tail surfaces. 

The aerodynamic forces acting on a wing section depend on the displacement h/b 
of a reference point along the chord (measured in semichord lengths b), the angle of rota- 
tion a as well as on their first and second time derivatives, that is, li/b, a, h/b, and 
01 . (See sketch (a).) It can therefore be expected that the L.E. control rotation /3 and 
the T.E. control rotation 6 will be a function of these parameters; that is. 



where [C], [G], and [R] are square 2 by 2 matrices. The reference frequencies coj. 
and Wj. are introduced to maintain the nondimensionality of equation (1) and can be 
used as free gain parameters. The matrix [R] can be neglected since the aerodynamic 
inertia forces are relatively small at the usual range of flutter frequencies. Hence, the 
following basic control law will be used during the initial stages of this work (an extended 
control law is suggested later in this work) ; 
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( 2 ) 



The aerodynamic energy analysis shows (ref. 17) that the optimum values of [C] and 
[G] are insensitive to the value of o>p. Furthermore, the smaller cop is, the more 
effective the active controls become. Equation (2) can also be written as 



The values recommended in reference 17 for the matrices [C] and [G] have 
been adopted in the present work; that is, 

5.6 
-1.4 

1.5 

-0.4 0.1 



Changes are introduced in the values of Gj ^2 reduce the maximum L.E. control rota- 
tion, whenever necessary, and to bring it closer to the maximum T.E. control rotation. 

It is important to note here that imlike in reference 17, h/h is taken to be positive 
in the upward direction. This case requires a change of sign of the first column of each 
of the optimized matrices [C] and [G] which appear in reference 17. Modifications 
of the basic control law are introduced within the section "Supplementary Investigation" 
presented later. 


APPROACH 

Two methods are common in gust analysis - the discrete gust approach and the 
continuous gust approach. There are some distinct advantages to the continuous gust 
approach. It is not biased toward a specific gust frequency or gust shape since it uses 
atmospheric turbulence information as input. Furthermore, oscillatory aerod 3 mamic 
coefficients (including control-surface coefficients) are used and these are readily avail- 
able for both compressible and incompressible flows. Its main disadvantage lies, how- 
ever, in the fact that statistical quantities are involved as response output. These 
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quantities do not lend themselves to an easy understanding of the physical side of the 
problem. The discrete gust approach enables one to follow the response of the aircraft 
in the time domain and thus promotes the understanding of the physical aspects of the 
problem. However, it has its own disadvantages. It has a distinct bias toward the gust 
input parameters, such as the gust frequency and gust velocity (often stipulated by federal 
aviation authorities), and the a.erodynamic coefficients relating to control rotations or to 
unsteady compressible flow are not readily available. Nevertheless, the discrete gust 
approach has been adopted in this work since the physical aspects of the problem are 
considered to be very important. 


Mathematical Model 

The whole aircraft is treated as one dynamic system. It is assumed that the motion 
of the aircraft can be described by its rigid-body translation, rigid-body rotation, and the 
superposition of up to 10 elastic modes. The wing and horizontal tail are each permitted 
to accommodate up to 10 activated strips. Time lag in gust encounter between the wing 
and the horizontal tail is incorporated in the dynamic system. 

Aerod 3 nnamics 

The nonsteady aerodynamic forces arising due to a step deflection of either the L.E. 
or T.E. control surfaces are derived in the appendix for two-dimensional flow and zero 
Mach number. Nonsteady strip theory aerodynamics (involving the classical Wagner and 
Kiissner functions) is employed throughout this work. The resulting aerodynamic forces 
are modified by introducing approximate finite-span effects through the square-root fac- 
tor which yields the well-known wing-tip singularity. In addition, approximate compress- 
ibility effects are simulated through the introduction of the Prandtl-Glauert correction. 

Computer Program 

For any combination of activated control strips, the program provides the time his- 
tories of the motions of the airplane, elastic deformation of the structure, control -surface 
rotations, and the bending-moment loads. Flutter speed is determined from observations 
of the time history of the oscillation associated with the structural deformation. 

APPLICATION 
Description of the Aircraft 

The complete geometric and djntiamic data of two aircraft were made available to 
the writer for the purpose of the present investigation. These aircraft are the twin-boom 
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Arava STOL transport and the Westwind business jet transport. These two considerably 
different aircraft are briefly described in the following sections. 

Arava aircraft .- The Arava (fig. 1) is a light-loading, low-speed, low-altitude STOL 
twin-turboprop transport (maximum mass 6800 kg). The high wings are hinged (in bend- 
ing) at their roots and retained by struts which are attached to the wings at points lying 
ahead of the quarter-chord point around half the semispan. The twin-boom structure 
which carries the tail unit provides strong elastic coupling between wing and tail. In gen- 
eral, the aircraft is a highly elastic light transport. 

Westwind aircraft . - The Westwind (fig. 2) is the Israel Aircraft Industry modified 
version of the Rockwell Jet Commander and has heavier loading, high-speed and high- 
altitude capability. It is a twin- jet business transport with a cantilevered wing (maximum 
mass 9400 kg). The wing is clean with engines fitted at the rear of the fuselage and car- 
ries empty tip tanks (chosen configuration). 

Spanwise Location of Activated Strips 

The wing of each aircraft is diyided into 10 equally spaced strips as shown in fig- 
ure 3. Each strip can accommodate a pair of active controls (that is, L.E.-T.E. controls). 
A strip with active controls is referred to as an activated strip. The strips located along 
the horizontal tail are allowed spans equal to one-third and one-tenth of the horizontal-tail 
semispan of the Arava and Westwind aircraft, respectively. (See fig. 3.) All leading- and 
trailing-edge control surfaces are 20 percent chord. 

In an attempt to reduce labor, the effect of Wp on the gust response of each air- 
craft is investigated at only five chosen spanwise locations along the wing and one loca- 
tion along the tail. 


Gust and Flight Speeds 

The gust disturbance is assigned a 1 - cos shape and its peak value is set to 
9.1 m/s for the Arava and 15.2 m/s for the Westwind for flights at Vq (the design cruis- 
ing speed which is taken everywhere in this work as the flight speed, except for the flutter 
investigation). The flight speed V^. of the Arava and Westwind are 85 m/s and 253 m/s, 
respectively. 

Mention should be made here that the 1 - cos gust shape 3 delds more information 
regarding the response of the aircraft as compared with the step gust case previously 
employed (ref. 20) and is therefore believed to be advantageous. 

Calculations 

■f 

The time histories of the airplane and control system responses (that is, h, a, ft 
6, li, a, ft 6 , ii, a, and wing bending moment) and the flutter speed are determined 
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for a range of values of iOj^. Since the flight and gust velocities are kept constant (except 
for the determination of the flutter speed), the response of the control system is directly 
related to the value of Wp. In other words, the value of 6max> example, is directly 
related, under these conditions, to the value of ojp. 

The variation with cop of the peak root bending moment, for example, is equivalent 
and more meaningful when presented in terms of order to generalize the results 

obtained herein and permit their presentation in a condensed form, cross plotting of the 
response time histories is found to be necessary. The definition of some of the quantities 
involved is illustrated in the following sketches. (See sketch (b).) 



Sketch (b) 

To a first order, the change in responses at the first peak ^that is, (^^bimax’ 
Ahmax) is directly proportional to and also to the gust velocity Vg. These can be 
made independent of the gust velocity if consideration is made of quantities such as 

(^^b)max/(^b)max’ ^^max/hmax, and dcj^max- As already stated, the values of 6c 

and 6max are directly related to Wp and, in turn, can be made independent of the gust 


9 



velocity when represented in terms of 


or 


provided the same flight 


speed is maintained. 


Vg/Ve 



Additional details regarding the general representation of the cross -plotted results, 
for both flutter suppression and gust alleviation, are given together with the presentation 
of the results. 


Scope of the Anal3i;ical Investigation 

Almost all the work reported herein relates to the activation of a single strip at dif- 
ferent spanwise locations for each of the two aircraft. The results for each activated 
strip provide information regarding both gust alleviation and flutter suppression and are 
presented in a manner which enables some generalization to a multistrip activation (one 
example is included). Two control laws are investigated: the basic control law (pre- 
sented under "Preliminary Investigation") which utilizes the total response of the airplane 
and the extended control law (presented xmder "Supplementary Investigation") which pri- 
marily employs the elastic responses of the structure. It is hoped that by treating two 
dynamically and widely different aircraft, as is the case in this paper, results are obtained 
which might be applicable to a wide variety of aircraft. 

PRELIMINARY INVESTIGATION 

The preliminary investigation is aimed at gaging the problems associated with the 
activation of a single preoptimized strip at different locations along the semispan of the 
wing and horizontal tail of both aircraft. The magnitude of the control -surface rotations 
is varied by changing the value of the free gain parameter l/co^. The effects of these 
control rotations on the gust-alleviation characteristics and flutter speed of the aircraft 
is then determined. 

In order to gain insight into the effect of active controls on the gust response of each 
aircraft, a sweep is made of a single activated strip along the five chosen wing stations of 
each aircraft. However, in order to maintain a common basis for comparisons, the 
parameter cUj. at each station is assigned a value which yields a maximum absolute 
value of control deflection of 0.5 rad. A generalization which covers a range of control- 
surface deflections is subsequently made. 

Results and Discussion of Preliminary Investigation 

The results are presented and discussed in three groups. These groups relate to 
the effect of active controls on the maximum bending moment of the wing, on the fuselage 
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accelerations, and on the flutter speed. Each such group contains the results available 
for both the Arava and the Westwind transports. 

Bending-moment effects .- Figures 4 and 5 show the effects of the spanwise location 
of the activated strip on the maximum bending moment of the wing for maximum control 
deflections of 0.5 rad. Mention should be made here that the maximum wing bending 
moment for the Arava transport is located at the inboard boundary of strip 5 where the 
struts join the wing. This location is unusual for conventionally built aircraft which have 
cantilevered wings and therefore develop maximum bending moment at the root of their 
wings. The Westwind, which is of conventional design, shows maximum values of bending 
moment at the inboard boundary of strip 10. 

From figures 4 and 5 the following points emerge: 

(a) The active strip, when located at the tip region of the wing is effective in reduc- 

ing the positive bending moment but, at the same time, leads to an increase in the negative 
values of the bending moment. This may be explained in part by the fact that the external 
forces (fig. 6), which are comprised of aerodynamic forces F^ and inertia forces 

Fj (that is, reverse mass acceleration forces), act through an inboard section of the wing, 
whereas the control-surface reaction Fc to these forces acts through the activated strip 
located at the tip region. In addition, the control force F^ reduces the acceleration of 
the wing in proportion to AFj. Such an arrangement promotes the bending of the wing 
and may lead, eventually, to a large bending moment. (Note that the 1 - cos excitation 
is allowed for the duration of a single cycle only.) Figure 7 illustrates the case where 
the activated strip is located at the root of the wing. It can be seen that the contribution 
of F(. to the bending moment is negligible (in this case) whereas the reduction of AFj 
leads to an increase in the maximum values of the positive bending moment. The best 
location of F^ to relieve the bending moment (for equal magnitudes of alleviated positive 
and negative peaks) is at the region of application of Fj^g (that is, at the external forces 
center). When F^ is large and is brought about by a single activated strip located at the 
external forces center, it might shift the maximum bending-moment section toward the 
outboard of the wing (since the external forces are distributed over the whole of the wing). 
This case, however, is not realistic since the control deflections required to produce such 
magnitudes of F^ are beyond the practical range of a single strip. Figures 4 and 5 show 
that the external forces center of the Arava is located at around 55 percent of the span 
(just inboard of strip 5) whereas the external forces center of the Westwind is aroimd the 
35-percent span region (just inboard of strip 7). The Arava case is somewhat more com- 
plex when compared with the Westwind because of the influence of the strut's reaction. 
Nevertheless, the location of its external forces center can be determined from figure 4. 

(b) The effectiveness of the activated controls is very different in the two aircraft 
treated; the Westwind shows a large reduction in the positive bending moment whereas the 
Arava shows a smaller reduction. The reason for this difference is threefold: 
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(1) The strut introduces a basic difference which pushes the external forces center 
to a region very near the maximum bending-moment section and hence reduces the effec- 
tiveness of the activated control for positive bending moment. 


(2) The bending moment is a linear function of the gust angle Vgy'V when the same 

flight speed is maintained. The activated control force is essentially a function of 6c 
which designates the control deflection at the time the bending moment of the unactivated 
aircraft reaches its first positive peak value ((^b)max)' effect of the L.E. control 

on the maximum bending moment is negligible compared with the effect of the T.E. con- 
trol. Hence, for each station, the reduction in the positive peak value of the bending 


moment 

(^“b)„ax/(“bLax) 

that large values of 


as a fraction of the positive peak bending moment (that is. 


will be approximately proportional to 


6. 


g/'^c 


Vo./V, 


. This relation implies 


Vg/Vc 


should lead to larger reductions in the positive peak bending 


moment. Hence, the effectiveness of a given 6g value depends on the value of Vg^ Vc- 

The smaller VgyVc is, the more effective the activated control becomes for a given 5q. 
The value of Vgy^Vc for the Westwind is around 0.06 whereas the corresponding value for 
the Arava is around 0.11; thus, the increased effectiveness of the Westwind is explained. 

In general, faster aircraft have smaller Vg^Vg values and this gives them an advantage 
over slower aircraft. 


(3) The constraint on the control-surface rotation applies to the maximum value of 
rotation. It therefore follows that the larger 6c/'5inax is, the more effective the active 


control becomes 


-^max 


-’max 


in all cases being maintained through appropriate 


variation of Gj2j‘ The values of 6cy/6max for the Westwind are around 0.95 whereas 
the corresponding values for the Arava are aroxmd 0.65. This phenomenon further 
explains the reduced effectiveness of the Arava control as compared with the Westwind. 


A logical outcome of this discussion is the representation of 


as a function of 


Vg/V, 


and 6c j 6jnax a function of 


Vc/V, 


max/\''‘J^max 

at the different stations 


along the span of the wing. This representation is meaningful provided the activated con- 
trols do not act to promote increased bending response at later time intervals, that is, 
when the control forces act to the outboard of the external forces center. Hence, the 

results pertinent to (^^b)j^axy/(^'^)max ^ function of ^ at strip 1, for exaihple 


12 



I 


(that is, at the wing tip), are meaningful provided the control surface at strip 1 is activated 
in conjvinction with other active strips located along the span of the wing in a manner which 
insures that the resultant control force acts through the external forces center. Figures 8 

and 9 show this variation of (aM^) /(^h) ^cAmax as a fmiction of ? — 

for both the Arava and the Westwind transports. The bending-moment curves show the 
expected linear relationship between (^^b)max/(^^^max y Jy (except for very 

large bending- moment reductions which shift the activated maximum bending-moment 

/ ®c 

curve toward the origin) . The curves of 5{./6jyj3x against — for both aircraft 

^g/^c . 

show some slight nonlinear effects. 

As already stated, the parameter ScAmax i® very important since it has a strong 

effect on the effectiveness of the control rotations. It seems reasonable to expect that 
6c/6max will, in general, be a function of both the gust frequency and the natural fre- 
quencies of the aircraft. The gust frequency of the Westwind is around 5 Hz and is very 
near the first wing bending frequency (==6.5 Hz). In the Arava case, the gust frequency is 
around 1.7 Hz whereas the first wing bending frequency is arovind 6.5 Hz. 

The effect of the values of airplane frequencies on the value of 6c Amax illus- 
trated in the following discussion. Let ^c/^xa 2 ix. be equal to unity. This means that 6 
varies in phase with the bending moment. However, the bending moment itself varies 
almost in phase with the external forces acting on the aircraft. It therefore follows that 
by attempting to obtain 6c/6max = 6 will vary in phase with the external forces. 

Equation (2)- shows that if the extreme cases where Wf is large are excluded, the second 
term in this equation can then be made dominant. This result implies that 6 varies in 
phase with the velocity of the response. Hence, the stated object of having 6 vary in 
phase with the external forces is equivalent to having the response velocity vary in phase 
with the external forces. Such a relationship exists only when the external forces excite 
the structure at one of its natural frequencies. The frequency of the rigid- body transla- 
tion is zero; thus, the first wing bending frequency is left as the lowest natural frequency 
which might substantially affect 

At this stage it is appropriate to draw attention to the fact that some rigid-body 
excitation will take place even when the gust frequency is identical to the wing natural 
frequency. This is particularly true since the gust excitation forces used in the present 
work (that is, 1 - cos) do not have a zero mean value. However, the large elastic 
response, in the latter case, is dominant and subsequently leads to large 
ues. Mention should also be made that the rotation of the controls introduce external 
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forces which lag both the control rotation' and the response velocity. Therefore, these 
considerations relate to general trends only and permit some deviations whenever the 
secondary effects become larger. 

Drawing attention to the gust frequency, it should be remembered that the Federal 
Aviation Regulations (F.A.R.) base this frequency on a 25-chord wavelength. Therefore, 
the faster the aircraft, the higher the value of the gust frequency. These higher gust fre- 
quencies, acting on faster flsdng aircraft, are thus brought closer to the first wing bending 
frequency and are generally responsible for the development of large gust loads. It is 
therefore interesting to note that the bending- moment alleviation using active controls is 
most effective ^through'the values of Sc/Smax 3-nd when the gust loads are 

largest. Similar conclusions can be obtained by using the continuous atmospheric tur- 
bulence approach since most of the turbulence frequencies lie within a frequency region 
which is strongly dependent on the flight speed. 

The curves presented in figures 8 and 9 can readily be used to yield results for dif- 
ferent values of 6jnax ^.nd Vg. Take figure 9 (Westwind) as an example and determine 
6c so that 5max exceed 1 rad at strip 1 when Vgy/Vc = 0.2. Note that the slope 

of a straight line (through the origin) relative to the ordinate is given by . In this 

^ Vg/Vc 

example = 5 and the straight line shown in figure 9 crosses the curve at point A 

Vg/Vc 

to give 7 — = -4.4. With this value of — ^ can be deter- 

Vg/Vc Vg/Vc ^ °^max/v o^max 

mined from the same figure, which yields the value of -0.2 (that is, 20 percent reduction). 
The advantages of the graphical representation are now clearly evident. 

Fuselage acceleration e ffects.- The various fuselage acceleration effects are pre- 
sented in this section. 

Center of gravity: The effect of the spanwise location of the activated strip on the 
center-of-gravity acceleration, for 6max = 0.5 rad is shown in figures 10 and 11. Here 
again it can be seen that the activated strips at the wing-tip region are relatively ineffec- 
tive in reducing the center-of-gravity accelerations. The reason for this ineffectiveness 
originates from the fact that the activated strip at the tip region reduces the local wing 
acceleration and is unable to reduce the center-of-gravity acceleration because of the 
promotion of wing bending, as explained in an earlier section of this paper. The best 
spanwise location of the activated strip for maximum reduction in the center-of-gravity 
acceleration can be seen to lie at the wing- root section (strip 10), where the local aerody- 
namic coefficients per unit span are largest (because of the finite-span effect) and where 
direct bending effects are negligible. 
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In a manner similar to the bending-moment case, one can define a new value of 
6{. which corresponds to the control-surface deflection when the center-of-gravity 
acceleration of the vmactivated aircraft reaches its first maximum and can represent 

(Ahcff) /(hcff) as a function of — ^ — and as a function of ^ — . 

^ ‘^S^max/^ ‘^SVax ^g/^c ' 

(See figs. 12 and 13.) The points previously discussed of the differences in effectiveness 

^c / 

between the Westwind and the Arava transports due to the different n — and 6c/®max 

Vgy'Vc 

values apply also to the center-of-gravity accelerations and confirm once again the 
advantages of fast aircraft over slower ones in regard to gust-alleviation problems. 

It might be of interest to note that the peak value of h^g is almost in phase with 
the peak value of the bending moment; thus, almost identical 6(> values are obtained. 

/ ^c 

For the Arava case, the curve representing /6niax ^ function of - for the 

center-of-gravity accelerations is, in fact, identical to the one obtained from the bending- 
moment variation and is therefore not included in figure 12. 

Acceleration at other stations along the fuselage: It is important to reduce the 
levels of accelerations not only at the center of gravity of the aircraft but also at differ- 
ent stations along the fuselage. Three points are chosen to indicate the accelerations of 
the fuselage. One of these points coincides with the center of gravity of each of the air- 
craft and has already been discussed. A second point, designated as point 1, is located 
along the axis of the fuselage aft of the center of gravity, and a third point, designated as 
point 2, is located along the same axis, forward of the center of gravity. For the Arava 
transport, points 1 and 2 are located at 3 m aft of and 3 m forward of the center of gravity. 
For the Westwind transport, points 1 and 2 are located 2.58 m aft of and 6.73 m forward 
of the center of gravity, respectively. Since both of the aircrafts vinder investigation are 
transport planes, a reduction in the accelerations at these three points, which span the 
length of the fuselage, has a strong influence on the ride comfort of the aircraft. 

Figures 14 and 15 show the accelerations at points 1 and 2 for the two aircraft for 
maximum control rotations of 0.5 rad. The similarity between the accelerations at these 
points and the center-of-gravity acceleration is apparent and it shows the increased effec- 
tiveness of the inboard stations over the outboard ones. It can be seen that the accelera- 
tions at point 2 (forward of the center of gravity) are smaller than the center-of-gravity 
accelerations whereas the accelerations at point 1 (aft of the center of gravity) are larger 
than those of the center of gravity. Thus, the existence of a pitch-down acceleration is 
indicated. 
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Pitch acceleration; Figures 16 and 17 show the variation of the pitch-down accel- 
erations around the center of gravity (rigid-body acceleration). It can be seen that the 
activated wing strips have little effect on the pitch accelerations. 


Flutter speed e ffects.- It is useful to relate the flutter speed results with the gust 
response results through the representation of AVf /Vf as a function of j — , where 


AVf represents the increase in the flutter speed due to the activated control, Vf repre- 
sents the flutter speed of the unactivated airplane, and Sq relates to the value of 6 
where the bending moment is maximum. For flutter speed calculations, the value of 
is irrelevant. However, it serves as a measure of and can be used to define the 

/ ®c 

value of Sjnax through the use of 6c/^max ^ function of — . 

^g/^c 


The flutter speeds for the unactivated Arava and Westwind aircraft are 233 m/s and 
320 m/s, respectively. The high flutter speed of the Westwind leaves a very small margin 
for AVf before crossing into the supersonic regime. Therefore, some changes in con- 
figuration were made (by essentially increasing the mass of the wings) which reduced the 
flutter speed to 240 m/s. Hence, the Westwind flutter results relate (unfortimately) to a 
different configuration and therefore, special computer runs were made to determine the 
/ 

variation of 6c/6jjjax 7 — changed configuration. It should further be 

Vg/Vc 


mentioned here that the activation of the control surfaces at high speed introduces, in the 
very low frequency range, an instability which looks like dynamic divergence. To counter- 
act this instability, Gjj is allowed a value of -0.4 for the Arava transport and a value of 
-1 for the Westwind. Figures 18 and 19 relate therefore to somewhat different control 


law parameters. Figures 18 and 19 show the variation of AVf /Vf with 


6, 


Vg/Vc 


and 


^c/^max 


with 


6. 


Vg/Ve 


with the spanwise location of the activated strip. It can be seen 


that, as expected, the most effective location of the activated strip for flutter suppression 
is at the tip region. For purposes of illustration, imagine that maximum control- surface 
rotations of 0.18 rad at Vg are allowed. For the Westwind, this assumption yields a 


value of 


^max 

Vg/Vc 


aroimd -3. Figure 19 (strip 1) shows that for a slope of -3 with respect 


to the ordinate, a value of about -3 is obtained for 


Vo-/V, 




With this value for 


Vo-/V, 


figure 19 shows that an increase of about 70 percent in flutter speed is obtained. 
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Mention should be made here that significant values of AVf^Vf bring the flutter 
speeds into the supersonic regime while subsonic aerodynamics is used. However, 

AVf/Vf should be considered in those cases only as an indication of beneficial effects 
and should not be used quantitatively. It is very interesting to note that the introduction 
of Gji leads to large increases in the values of the flutter speed for both the Arava and 
Westwind transports. Further increase in the Gu values leads to a further increase 
in the flutter speeds. The parameter G^ is assigned the smallest value which stabi- 
lizes the divergence instability to the extent that the flutter speed can be determined. The 
effect of the parameter Gn, which is associated with the L.E. control, stresses once 
again the importance of the L.E. control. 


Application of resu lts to multistrip activation .- It can be seen that a single activated 
strip cannot be effective for flutter suppression without leading to a degradation in the 
bending- moment values. Furthermore, two or more activated strips are required to 
tackle both the flutter- suppression and gust-alleviation problems (that is, bending moment 
and center-of-gravity accelerations). These activated strips should be placed at the wing 
tip (for flutter suppression) and in the wing-root regions (for center-of-gravity accelera- 
tion) in a manner which insures that the resultant control force passes through the exter- 
nal forces center (for bending- moment alleviation). The representation of the results as 
Sf. 

a function of ^ — can also be used to estimate the joint effects of two or three neigh- 

Vg/Vc 

boring activated strips. This case can best be illustrated by the following example: 


Assume that it is desirable to find the effect of three activated strips located at 
strips 4, 5, and 6 of the Arava wing, maximum control deflections of 0.5 rad being allowed 
at each station. Since the resultant force of the three activated strips will pass aroxind 
strip 5, it will be equivalent to a single force caused by a single activated strip located at 
strip 5, having a maximum control deflection of 1.5 rads (that is, 0.5 x 3). This condition 

leads to = 13.6 which yields the value of 7 — = 9.9. (This value is the same 

Vg/Ve Vg/Ve 

for both the bending moment and fuselage accelerations. See fig. 8.) The gust-alleviation 
and flutter -suppression effects follow from figures 8, 12, and 18 by using this value of 
5c 

— . These figures jdeld the following values: 

Vg^Vc 

(am^) 

— — = -0.3 (that is, 30 percent alleviation) 

\^t))max 
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— — ■ = -0.43 (that is, 43 percent alleviation) 

The ratio AVfy/Vf was not calculated. Figure 20 shows the results obtained from a spe- 
cial computer run for the preceding three simultaneously activated controls. The agree- 
ment is seen to be extremely good; thus, the suggested procedure is confirmed. 

Activation of a co ntrol strip along the horizontal tail.- Initial attempts to place an 
activated strip on the horizontal tail have led to a severe static longitudinal instability. 
This condition indicates that some changes will have to be introduced into the control law 
to permit the location of an activated strip (with no adverse effects) on the horizontal tail. 
This point is discussed in a subsequent section. 

Summary of the Preliminary Investigation 
The main points emerging from the preliminary investigation are as follows: 

(1) A single activated strip promotes the bending of the wing when located at either 
the outboard region of the wing or at its root region. 

(2) The optimum location of the activated strip for bending- moment alleviation is 
around the external forces center region. 

(3) An activated strip located at the inboard region of the wing is most effective in 
reducing the center -of- gravity acceleration. 

(4) An activated strip located at the outboard region of the wing is most effective for 
flutter suppression. 

(5) The value of 5(,/6max tends to unity when the gust frequency approaches one of 
the wing's natural frequencies. 

(6) The representation of (^^cg)^ /(hcg)^^, and AVf/Vf, 

/ 5c 

6c/6jjiax as a function of — — allowed the generalization of the results over a wide 

Vg/Vc 

range of conditions. 

SUPPLEMENTARY INVESTIGATION 

In the following section, the problems arising from the bending- moment degradation 
when locating a single activated strip at the wing-tip region are investigated and attempts 
are made to locate an activated strip on the horizontal tail with no adverse effects on the 
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longitudinal stability of the airplane. These attempts lead to the formulation of an 
extended control law which is subsequently applied to the basic problem of flutter 
suppression and gust alleviation. 


Extended Control Law 

It has already been mentioned, in previous sections, that the location of an activated 
control strip at the tip region of the wing leads to an increase in the absolute bending- 
moment values. This increase is clearly a result of the coupling (through the activated 
aerodynamic forces) between the rigid-body movements, due to the gust forces, and the 
wing elastic modes. These coupling effects promote the bending of the wing and require 
the addition of one or more active strips at the root region of the wing to force the resul- 
tant control force to pass through the external forces center. Furthermore, the rigid- 
body motion and these coupling effects introduce large wing responses due to gust, pro- 
duce large control deflections, and thus limit Wj. to a relatively large value. This 
limitation of Wj. may eventually lead to a degraded performance of the activated strip 
at the flutter frequency since, in general, the effectiveness of the activated control 
increases with 1 /wj. (except for those cases where Wr is very small, as indicated in 
ref. 17). Hence, a control law which can either eliminate or reduce the rigid-body input 
to the control law or introduce some controlling elements to this latter input might prove 
to be superior. Such a reduction can be achieved by the following control law: 



The values of h and a at the reference section are denoted in equation (5) as h^. 
and 

The first two elements in equation (4) control the relative motion of the wing where- 
as the last two elements control the reference section movement. By using equation (5), 
equation (4) reduces to 
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In order to reduce the rigid-body contribution to the control- surface rotations, the root 
section is chosen as the reference section for control systems located along the wing. 

It should be observed that when 

equation (6) reduces to the basic control law (eq. (2)). This result implies that, in gen- 
eral, the reference section of an elastic structure can be controlled to a limited extent. 

If too large forces are applied on an, elastic structure to control the motion of the refer- 
ence section, dynamic instability might result. If, however, the magnitude of the forces 
which control the motion of the reference section is limited to values smaller than those 
defined by o>p = 2 > the stability of the system is expected to be maintained. It will 

therefore be postulated that stability will generally be maintained provided 

^r,2 = ‘^r (7) 

Hundreds of computer runs were made for the two aircraft by using a variety of 
combinations between ooj. and 2 while maintaining equation (7). During those runs, 
no instability developed; thus, some measure of confirmation is provided to this postula- 
tion. The results obtained by using this extended control law are presented in the follow- 
ing sections. 


Results and Discussion of Supplementary Investigation 

The results obtained through the use of the extended control law are presented and 
discussed in three groups in much the same way as those relating to the ordinary control 
law. Each such group contains the results obtained for both the Arava and the Westwind 
transports and is accompanied by a discussion of these results. In all cases the value of 
2 assigned large values compared with (at least one order of magnitude larger) 

in order to determine the effect of the extended control law near the extreme end of its 
range. 

Bending-moment effects .- Figures 21 and 22 show the effects of the spanwise loca- 
tion of the activated strip on the maximum bending moment of the wing that are appropriate 
for a maximum control deflection of 0.5 rad. It should be observed that the results relat- 
ing to the root section of both aircraft (strip 10) are not presented. The reason for this 
omission originates from the fact that since the root section of the wing has been chosen 
as the reference section, the extended control law reduces, in this particular case, to the 
ordinary control law which has already been treated. It should also be noted that the 
results relating to the activation of strip 8 of the Arava transport are not presented 
because of the extremely small values of required to produce moderate control 


20 



rotations. The results relating to strip 7 of the Westwind are very similar to those of 
strip 6 and have not been included here. The following discussion therefore relates to 
strips 1, 3, and 5 of the Arava and to strips 1, 4, and 6 of the Westwind. Figures 21 
and 22 indicate that by eliminating the rigid-body velocities from the control law, the 
activated strip can be located at the tip region of the wing without promoting bending. 

The optimum location of the strip appears to be somewhat inboard of the wing-tip strip 
(strip 1). This result can clearly be seen in figure 22 which relates to the Westwind, 
where the results for strip 4 seem to indicate the optimality of this strip. Figure 21 
which relates to the Arava, however, yields almost identical reductions in bending moment 
for both strips 1 and 3 and indicates that the optimum lies arovind strip 2. The reason 
for the inboard location of the optimum station is due to the finite -wing effects (elliptic 
distribution at the tip) and to the tapering of the wing (Westwind only) toward its tip. 

A comparison between the results relating to the ordinary control law (figs. 4 and 5) 
and those relating to the extended control law (figs. 21 and 22) with a maximum control 
deflection of 0.5 rad in both cases shows the superiority of the extended control law. As 
an example, strip 4 of the Westwind shows almost a total reduction of the bending moment 
when using the extended control law, whereas the same strip with the same maximum con- 
trol deflection yields only a reduction of about 60 percent. It should also be recalled that 
the 60-percent reduction in bending moment can be obtained provided other activated strips 
are located aroimd the root section of the wing to insure that the resultant control force 
acts through the external forces center. It should further be recalled that the activated 
strips located inboard of the external forces center introduce detrimental effects regard- 
ing this positive bending- moment reduction and thus yield an overall bending-moment 
reduction which is smaller than 60 percent. It can therefore be seen that the basic con- 
trol law (which enables the reduction of the bending moment, in this case, through the 
activation of more than one strip) is distinctly inferior to the extended control law in 
regard to bending-moment effects. 


If a strict comparison on the basis of a single strip only is made, then the results 
relating to station 4 by use of the extended control law should be compared with those 
relating to station 7 (which almost coincide with the external forces center) by use of the 
basic control law. Figure 5 shows that for a maximum control deflection of 0.5 rad, the 
bending-moment reduction of the Westwind is around 35 percent; thus, the effectiveness 
of the extended control law is stressed. The extended control law is now investigated 
over a range of maximum control deflections and the results are summarized in figures 23 
and 24 by using the same representation as adopted for the ordinary control law, that is. 


^c/^max as a function of 


Vg/V. 


and as a function of 




max 


g; 


:/Vb' 


These remarks regarding the effectiveness of the extended control law can be seen to be 
appropriate over the entire range of control deflections. At the lower range of these 
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deflections, the relative effectiveness of the extended control law is somewhat diminished 
because of the reduction in the values of 6c/^max- reason for this reduction ^in 

6c/5max) relative values of the various terms in equation (6). As already stated, 

assigned large values compared with Wj.. As a result, the third term in equa- 
tion (6) is very small in most cases. The first term in equation (6) has not been modified 
and is identical to the first term in equation (2). Thus, the second term of equation (6) 
represents the main difference between the basic control law and the extended control law. 
At this stage it becomes apparent that the second term in equation (6) must be substantially 
smaller than its equivalent counterpart in equation (2) (for identical values of o)j. and 
[G]) because of the smaller values of the relative velocities ii and a (employed in 
eq. (6)) compared with the space velocities h and a (which are employed in eq. (2)). 
Thus, in order to obtain identical values of 6 through the use of equations (2) and (6), 
the values of coj. must be substantially reduced in equation (6) in order to compensate 
for the reduced response velocities ii and a. If it is remembered, however, that the 
second term gives rise to aerodynamic damping forces, it follows that the stated reduc- 
tion in coj. leads to larger values of the aerodynamic damping coefficients [which are 


proportional to — [G]j. It can therefore be stated that at the lower range of — ~ — , 

^r / ^g/^c 

the values of Wj, are not small enough to render the second term in equation (6) domi- 

®c 

nant. At the upper range of — , this second term increases and leads to control 

Vg^Vc 

rotations which are almost in phase with the relative response velocities ii and a. For 
the remainder of the discussion, the argument regarding the conditions which yield large 
values of 5c/^max follows along lines identical to those discussed earlier. 


Figure 23 which relates to the Arava transport with the extended control law shows 
-'max 


larger values of 


6c/6r 


at the upper range of 


Va- Vr 


compared with those obtained 


through the use of the basic control law'(fig. 8). The reason for this increased effective- 
ness in 5(./6jnax extended control law yields large damping coef- 

ficients which introduce large damping forces directly into the elastic modes and almost 
completely ignores the rigid-body velocities. These large damping forces lead to changes 
in the phase relationship between the exciting force and the response of the wing in such a 
direction that the response tends toward 90° phase lag, irrespective of the frequency of 
excitation. Such a phase lag is beneficial for the increase in 6c/6maxj discussed ear- 
lier, and leads to the rapid decay of the response of the elastic modes. 
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It is interesting to note that it is almost impossible to increase 


— (by reduc- 

Vg/V, ^ 


ing beyond the value which almost completely reduces the bending moment. This 
condition indicates that the activated control strip has introduced very large damping 
forces. Furthermore, within the range where increases, there is a nonlinear 

relationship between the activated control effectiveness and the maximum control deflec- 
tion. Thus, if the maximum control deflection is allowed to double its value, the effec- 
tiveness of the activated strip in reducing the bending moment will increase by a factor 
larger than 2. To illustrate this point, consider point A in figure 24 (relating to strip 1). 

5 6 

This point corresponds to = -4.88 and — 7 — = -2. The point designated as 

Vg/Vc Vg/Ve 


point B corresponds to = -9.76 and 


/Vc 


Vg/Ve 


= -9.3. In passing from point A to 


point B, 


^max 

Vg/Ve 


is allowed to double its value; in so doing, 


Vg/Ve 


increases by a fac- 


tor of 4.65. The resulting alleviation in the bending moment can be seen to increase from 
23 percent (at point A) to 61 percent (at point B). This nonlinear dependence. on 6 max 
gives, once again, an advantage to high-speed aircraft because of their relatively large 


values of 


v„/v, 


It should be stated however that assumes small values when yielding interme- 
diate control-surface rotations leading to large values of (eq. (3)). Reference 17 

indicates that a deterioration in stability might develop at very large values of 
Subsequent work (ref. 21), however, indicates that some simple modifications can be 
introduced into the control law which may maintain the absolute stability of the control 
law even at very high values of 


Fuselage acceleration effects .- Figures 25 and 26 show the effects of the spanwise 
location of the activated strip on the center-of-gravity acceleration for maximum control 
deflections of 0.5 rad. It can be seen, once again, that the optimum location of the acti- 
vated strip for the maximum reduction of the center-of-gravity acceleration is near the 
root region of the wing. It should be remembered that at the root region itself, the 
extended control law reduces to the ordinary control law. In addition, at wing sections 
very near the root section, the values of cOp required to produce reasonable control 
rotations become excessively small. Comparisons between figures 25 and 26 and fig- 
ures 10 and 11 indicate that, in general, only minor differences can be detected between 
the two control laws. In those cases where the extended control law yields larger values 
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of ^c/^max> above-mentioned differences are favorable toward the extended control 
law (on the basis of limited control deflections). These remarks are mainly applicable 
to the Arava transport which shows relatively small values of ^c/^max when the basic 
control law is used. Figures 27 and 28 confirm these stated trends in the representation 

of (aKpo.) /fhpD-') as a fvmction of — ^ — and 5p/6mav as a function of — 7 — 

V ‘^S^max/V ‘^g^max Vg/Vc ' Vg/Vc 

/ ^c 

(for both aircrafts). It should be noted that the variation of 6c/5max with — for 

Vg/Vc 

the Arava transport is identical to the one shown in figure 23. Figures 29 and 30 repre- 
sent fuselage accelerations at points 1 and 2 that are appropriate for maximum control 
deflections of 0,5 rad. 

Flutter speed effects.- By following identical lines as in the case of the basic con- 

/ ^c 

trol law, the representation of AVf /Vf as a function of ■. — is shown in figures 31 

Vg/Vc 

and 32. For comparison purposes, Gjj is assigned the same values in the supplemen- 
tary investigation as those used in the preliminary investigation. The large advantages 
of the extended control law can immediately be seen upon comparison with the basic con- 
trol law results. (See figs. 18 and 19.) This comparison indicates that there is a very 
close relationship between the flutter- suppression control law parameters and those 
relating to bending- moment alleviation. It is observed that the activated strips are most 
effective in suppressing flutter, when located at the tip region of the wing, at strips which 
are outboard of the optimum locations for maximum bending- moment reduction. It is 
observed (figs. 31 and 32) that at the wing tip, flutter is suppressed with very small 
control -surface rotations of a single activated strip. The Arava transport requires 
6 max ~0.07 rad whereas the Westwind transport requires S^ax ~0.2 rad. The higher 
gains required for substantial bending-moment alleviation will also result in increased 
flutter-suppression effectiveness. 

It should be noted that most of the wind-tunnel tests which are performed with the 
purpose of investigating the effects of active controls on the flutter speed relate to models 
attached to a rigid support. (See refs. 11, 12, 14, and 15.) This condition is equivalent to 
using the extended control law, as presented in this paper. Hence, large differences may 
be observed between "free flying" models and the "clamped" ones. 

Finally, it is important to note that the use of the extended control law either elimi- 
nates the dynamic divergence instability or else considerably reduces its effect. It seems 
therefore that this dynamic divergence instability arises because of the coupling between 
rigid-body and elastic modes. Best results are obtained when a>j. 2 is made very large 
and, as a result, flutter-suppression effectiveness is increased and control rotations are 
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reduced. A reduced effectiveness can be observed in figure 31 (strip 5). Here ^^,2 = 2.5 
whereas in all other cases oJj. 2 assumes larger values (either 20 or 100). This value of 

6c 

CO- o explains the relatively large values of — appearing in figure 31. It appears 

\j'/c 

therefore that the complete elimination of the contributions of the rigid-body motions from 
equation (4), by leaving only the first two terms, is likely to yield best results for both 
flutter suppression and the elimination of the dynamic divergent instability; that is, 



Activation of a contro l strip on the horizontal tail . - It has been pointed out (ref. 17) 
that the two-dimensional control law, derived through the use of the aerodynamic energy 
concept, activates the controls in such a way as to counteract any lift build-up on the strip 
and provide large damping forces. This reduction in lift build-up eventually leads to both 
a reduction in the bending moment acting on the wing and a reduction in the acceleration 
sensed at the center of gravity of the aircraft. However, when placing the activated strip 
on the horizontal tail of the aircraft, the resulting reduction in the lift build-up due to a 
rigid-body motion clearly reduces the effectiveness of the horizontal tail and leads to a 
deterioration in the longitudinal static stability of the aircraft. To counteract these detri- 
mental effects on the static stability, the following control law, which completely eliminates 
the responses of the reference section, is suggested by equation (8). 

For gust -alleviation work, the reference point is chosen along the wing-root section 
which corresponds to the center -of -gravity location of the airplane. This choice implies 
that a simple rigid-body pitch movement will lead to 5 = 0 and thus neutralize the acti- 
vation of the control due to this movement. The gains associated with h in the [C] 
matrix are zero whereas the gain associated with ii and a lead to dissipative forces 
which may introduce some sluggishness into the pitch movement of the aircraft with no 
subsequent degradation in the static stability. For horizontal-tail flutter problems, the 
movement of the center of gravity of the aircraft is generally unimportant in controlling 
the instability and the activated strip should therefore maintain its effectiveness (since 
the elastic deformations are fully maintained in h. O', and their derivatives). 

The suggested control law is applied to both the Arava and the Westwind transports 
and the results relating to an activated strip located near the center-line region of the hor- 
izontal tail (strip 11) are shown in figures 33 to 36. Figures 33 and 34 show the fuselage 
accelerations and the maximum bending-moment variation following the activation of the 
tail strip with 5naax = 0.5 rad (the Arava transport results relate to 6rnax = 0.25 rad). 
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It can be seen that the pitch accelerations have been reduced (leading to a reduction in hx 
and an increase in h2), with particular reductions in the absolute peak negative values. 
Furthermore, a reduction can also be observed in the center-of-gravity acceleration and, 
here again, the largest reductions relate to the peak negative values. The bending-moment 
variation shows a slight increase in the peak positive bending moment with a substantial 
reduction in the peak negative value of the bending- moment curve. Figures 35 and 36 
summarize the results relating to the variation of (^%)max/min/(^b)max/min’ 


(Ahcg)j. 

(AK2), 


cgymax/min/ (^cg)xnax/min 




'^max/min/(^l)max/min’ 

^c 

as fimctions of ; — for the two aircraft. 




max/ min / V^2ymax/ min 


Vg/Ve 


(The sub- 


script max/min denotes that either the maximum or minimum values for the ratio are 
taken.) 


No problems relating to horizontal -tail flutter could be simulated on the two air- 
craft treated and therefore, active flutter suppression of tail surfaces could not be con- 
sidered. However, when tail flutter problems exist, the tail root section is believed to 
provide the best location for the reference point. It can therefore be seen that the sug- 
gested control law provides the necessary means for flutter suppression and gust allevi- 
ation by using activated strips on the tail surfaces. 

Remarks on extend ed control law.- The supplementary investigation included the 
introduction of an extended control law, based on the elimination of some of the rigid- 
body contributions from the active control input. This elimination of the rigid-body 
movements cannot be regarded as the normal filtering process widespread in control 
theory. The Westwind case is best suited for the purpose of illustrating this point; both 
the gust frequency and the first natural wing bending frequency are approximately 5 Hz. 
Hence, if it is desired to filter out the rigid-body response to gust input, the first wing 
bending frequency must be filtered out as well. Thus the activated controls are prevented 
from reacting to the fluttering oscillation which takes place at ~5 Hz. On the other hand, 
the flutter frequency of the Arava is ~12 Hz whereas the bulk of the gust frequencies is 
=1.7 Hz. This value of the frequency allows the rigid-body responses to be filtered, and 
the free elastic modes are left to be controlled by the active system. In this latter case, 
however, one cannot directly tackle the gust -alleviation problem or take advantage of the 
flutter-suppression active system to control the gust loads (because of the filtering out of 
the rigid-body frequencies). It can therefore be seen that both the extended control law 
and the tail activation introduce new elements by considering the relative movement of the 
structure as inputs to the control system (as against the elimination of all the input sig- 
nals below a certain frequency). 

Finally, it should be mentioned that the two aircraft treated in the present work have 
negligible wing sweep angles. It is expected that when large sweep angles exist, rigid-body 



pitching angles may activate the wing-tip strips, based on the extended law, in much the 
same way as it activates the strips on the horizontal tail. It can thus be seen that some 
measure of rigid-body pitch control can be obtained from the activated strips located on 
the wing. A control over these pitching forces can readily be obtained by moving the ref- 
erence point backward along the root chord. If necessary, each wing strip can be allowed 
a different reference point (for the extended control law) along the root chord; thus, the 
complete elimination of the rigid responses is insured. 

Summary of Supplementary Investigation 

The main points emerging from the supplementary investigation are as follows: 

(1) The extended control law introduces large damping forces directly into the elastic 
modes while being only slightly affected by the rigid-body response of the aircraft. 

(2) The effectiveness of the activated strip using the extended control law is greatly 
increased for both bending-moment alleviation and flutter suppression. 

(3) The extended control law appears to be somewhat more effective than the basic 
control law for outboard wing sections in regard to fuselage accelerations. This point is 
true for outboard strip locations since it leads to comparatively larger values.of Sc/Smax* 

(4) The optimum strip location for maximum reduction in fuselage accelerations is 
at the inboard region of the wing (but clearly not on the reference strip). 

(5) The optimum strip location for maximum bending- moment reductions is in the 
tip region of the wing inboard of the tip strip. 

(6) The extended control law eliminates the need to define the external forces center. 

(7) The optimum strip location for maximum increase in flutter speed is at the wing- 
tip strip. 

(8) The extended control law, with the reference section defined as the tail root 
section, can be used to suppress tail flutter. 

(9) Activation of a strip on the horizontal tail with the wing chord taken as the ref- 
erence section can be used mainly to reduce the absolute negative values of the bending 
moment and center -of -gravity acceleration (if necessary). 

CONCLUDING REMARKS 

The effectiveness of activated leading-edge and trailing-edge control systems on 
flutter suppression and gust alleviation has been determined in the present work. Most 
of the results have already been summarized in previous sections of this work. However, 
it is felt appropriate to mention here some of the major points which emerge from the 
present work: 
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(1) The activated system introduces strong coupling effects between the rigid-body 
response of the aircraft (to gust inputs) and the wing bending modes. These coupling 
effects impose severe restrictions on the location of the activated systems along the span 
of the wing (their resultant should pass through the external forces center) when bending- 
moment loads are important. 

(2) The coupling between the rigid-body modes and elastic modes, introduced by the 
activated system, leads to large control deflections and to a greatly reduced overall effec- 
tiveness of the control system. This effect is particularly pronounced when the activated 
system is located at the outboard wing region, and it introduces severe limitations on 
bending-moment alleviation and flutter-suppression capabilities of the activated system. 

(3) The extended control law, which is based on the wing elastic deformation only, 
appears to overcome these limitations and leads to almost complete decoupling between 
the rigid-body responses, elastic responses, and the activated control forces. As a result, 
considerable improvements in both bending- moment alleviation and flutter suppression are 
obtained. 


(4) The extended control law permits the suppression of tail flutter without leading 
to losses in the static longitudinal stability of the aircraft. 


(5) It is found that the effectiveness of the activated control system ultimately 
depends on the value of the ratio between the maximum trailing-edge control deflection 

6max angle Vg/Vc, that is, . It is therefore very easy to present 


excellent gust-alleviation results for small 6max values whenever Vs/^c assigned 

' 6 . ' 


very small values to yield large values of 


^max 


V' 


Comparison between the performance 

6 , 


of various control systems should therefore be made only at comparable values of 


^max 

Vg/Vc' 


(6) It is found that activated control systems are more effective in gust alleviation 
and flutter suppression of fast flying aircraft when compared with slow flying aircraft. 

(7) Free fljdng wind-tunnel models might show a greatly reduced performance as 
compared with clamped models unless some form of an extended control law is used. 


(8) The basic control law is best suited for reducing the acceleration levels along 
the fuselage provided it is applied to activated strips located very near the root of the 
wing. Deterioration in bending- moment levels may result unless the resultant of these 
activated control forces is made to pass through the external forces center or a compen- 
sation based on the extended control law is made by placing an activated strip at the out- 
board end of the wing. 
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Many graphical results have been presented in the present work with the hope of 
contributing toward a better insight into the problem of flutter suppression and gust alle- 
viation and with the purpose of showing the strong interrelations that exist when attempt- 
ing to treat separately only some of the aspects associated with the wider problem. Addi- 
tional work is, however, required which incorporates the actuator dynamic characteristics, 
continuous gust inputs, three-dimensional aerodynamics, and performance at supersonic 
speeds. Nevertheless, it is very encouraging to discover the large effectiveness of the 
activated control systems in flutter suppression and gust alleviation. It also appears 
remarkable that a "utility" type, preoptimized control law, which is based on the concept 
of aerodynamic energy, brings about such stable and powerful control systems. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
April 27, 1976 
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APPENDIX 


DERIVATION OF THE GENERALIZED AERODYNAMIC FORCES OF A 
LEADING- EDGE AND TRAILING-EDGE SYSTEM DUE TO A 
UNIT STEP ROTATION OF THE CONTROL SURFACES 


In addition to the symbols given in the body of the paper, the following symbols are 
used in the derivation: 

e distance between midchord and leading edge of control surface ^also designated 

as fg) 

F(w) Fourier transform function defined by equation (A17) 

F(k) real part of Theodorsen fimction 

G(k) imaginary part of Theodorsen function 

h,h displacements defined in sketches (c) and (d), respectively 

I,Ij integrals defined by equations (A21) and (A24), respectively 

Jj integral defined by equation (A25) 

k reduced frequency, u)b/V 

aerodynamic forces defined in sketch (d) 


l,m,t;i 

P,Q,RJ 




p. Q. > constituents of L, M, T, P, Q, and R, respectively, comprised in j^ApJ 
3 ^ J ^ 3 I 


distances defined in sketch (d) 
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Qh/b»Qo!^ 

Q^>Q6 J 


Wh/b»Wo,;\ 

W^,W6 J 

OfJ,6 

(p{s) 

P 

Matrices; 

M 

[Ai],[A2] 

[bi].[b2] 

[D] 

[E] 

[Ho],[Hi] 

Notation: 

If 


generalized aerodynamic force along h/b, or, and 6, respectively 

ith generalized coordinate 
nondimensional distance, Vt/b 
jth Theodorsen T function 
time 

virtual work along coordinates h/b, o;, /3, and 6, respectively 

rotations defined in sketch (d) 

Wagner function 
fluid density 

aerodynamic coefficients matrix defined by equation (Allb) 
aerodynamic matrices defined by equations (A40) and (A41), respectively 

aerodynamic matrices defined by equations (A42) and (A43), respectively 

transformation matrix, defined by equation (Alla) 
aerodynamic matrix defined by equations (A44) 

matrices of aerodynamic coefficients defined by equations (A13) and (A14), 
respectively 

transposed matrix 
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Subscripts: 

s step-type variation 

6 parameter relates to 6 

Dots over S 3 unbols denote derivatives with respect to time. 

First consider the determination of the aerodynamic oscillatory forces acting on 
the leading- and trailing-edge system shown in sketch (c). Assume that the oscillatory 



Sketch (c) 

forces acting on a somewhat similar trailing edge and tab system (as used in refs. 22 
and 23) described in sketch (d) are known. Note the difference between the h, a, 
and 6 coordinates and the h, a, and 6 coordinates. The arrows in each sketch 
indicate the direction of positive displacements, forces, and distances. 


Undisturbed position 



The symbols L and M denote, respectively, the total lift and pitching moment 
and L is assumed to act through the quarter-chord point. The force acting on the 
trailing-edge and tab combination through the leading-edge point is denoted by P and 
the aileron tab moment is denoted by T. The force acting on the tab moment is Q. 
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The forces in sketch (d), their direction, and points of application are identical with those 
of Smilg, Wasserman, et al. (See refs. 22 and 23.) 

As a first stage of the following analysis, the oscillatory generalized forces acting 
along the h, a, |3, and 6 coordinates are determined through the application of the 
principle of virtual work. 

The virtual work Wj^ in the h coordinate is given by 

Wh = -b I L (Al) 

where b represents the semichord length. Denoting the generalized forces by sub- 
scripted Q's yields from equation (Al) 

6W. 

Qh /h = T- = -bL (A2a) 

5(h/b) 

Similarly, the virtual work W^, in the a coordinate is given by 


Wq, = Mo; - Lf jQ! 


or 


6Wq, 


The virtual work in the /3 coordinate is given by 

= -^M - T - + (L - P)£4/3 


and hence, 


^ = p(^2 - £4) + L«4 + T - M 


Similarly, 


6^ 


Wg = Q6 - R £36 


or 


(A2b) 


(A3) 


Qg = 


aw 

~65 


^ = Q - R£3 . . . 


(A4) 


33 


APPENDIX 


These equations can be condensed into the following matrix equation: 


r ^ 

%/b 


-b 

0 

0 

0 

0 

0 

Qq; 

< 

> = 


1 

0 

0 

0 

0 



f4 

-1 

1 

H-2 - ^4 

0 

0 

Qs 


0 

0 

0 

0 

1 

-^3 

-i 


fL^ 


M 

T 

P 

Q 




(A5) 


but the forces are of the form 


L 


b3 


h/b 

M 


b4 


a 

T 





< 

P 


b3 

[Ao]< 

z/b 

Q 


b4 


6 

R 

V J 




y/b 



L J 


(A6) 


where . |^AqJ is a 6 x 6 matrix and z and y define the hinge locations and are given by 

z = (f 2 - 
y = -^3^ 


Therefore, one can write the following equation based on the latter two equations: 
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h/b 


1 

0 

0 

0 

a 


0 

1 

0 

0 



0 

0 

1 

0 

^ z/b 

> = 

0 

0 

b 

0 

6 


0 

0 

0 

1 

y/b 


0 

0 

0 

col 



- 



J 


C— ~\ 
h/b 





L^J 


(A7) 


Comparison between the coordinates h, a, and 6 and h, a, /3, and 6 yields 
the following relations: 

h = -(h + £iQ! - 

5 = a? - j3 
^ = 13 

6 = 6 

These expressions reduce to 


r- ^ 

h/b 



-(l/b 

iijb 

0 


h/b 

a 

( 

) = 

0 

1 

-1 

0 

< 

a 

) 

/3 


0 

0 

1 

0 



6 

J 


0 

0 

0 

1 




(A8) 
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Substitution of equation (A8) into equation (A7) yields 


h/b 


-1 


^4 

0 




b 

b 



a 


0 

1 

-1 

0 


r 








h/b 



0 

0 

1 

0 

I 

a 


> = 



fo - Ha 



> 

z/b 


0 

0 

2 4 

b 

0 










6 



0 

0 

0 

1 


^ J 

y/b 


0 

0 

0 

col 
1 1 










(A9) 


Substituting equations (A6) and (A9) into equation (A5) yields the following relation between 
the generalized forces and the displacements of the L.E.-T..E. system 


where 


r 

^h/b 


r 

h/b 

Qflf 

) = irpoj2b4[D]'^|Ao][D]‘^ 

a 

< 

Q|3 

) 

/3 

Q 5 


6 

^ J 



(AlO) 


(Alla) 
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Letting 


H = 


^1 

Lg 


L4 

H 

H 

Ml 

M2 

Mg 

M4 

Mg 

Mq 

Tl 

T2 

T3 

T4 

T5 

Te 

Pi 

P2 

P3 

P4 

P5 

Pe 

Ql 

Q2 

Qg 

Q4 

Q5 

Qe 

Ri 

R2 

R3 

R4 

P5 

Re 


and using equation (AlO) yields 

= TTPCO 




where 


|Hol = 


[Ll] 



L, -A - Ml 


hi h \ ^1 


and 


N = 


Hi( 2 ,l) Hi( 2 , 2 ) 


where 


(Allb) 


(A12) 


(A13) 




^2 " ^41 


Hi(l,l) = -(Li ^ - L2 + Lg + L4 ^ 

Hi(l,2) = -Lg + Lg ^ 

Ha U - H. £. / e. 9-2 - 9. 

Hi( 2 ,l) = Ml ^ - M2 + Mg + M4 - ^(Li ^ - L2 + Lg + L4 — 

Hi(2,2)=M5-Mg^--^^L5-L6-^^ 


(A14) 
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Equation (A14) can be written as 

[»l] = [Ta][A][Tp] (A15) 

where 



In order to compute the forces due to control-surface step rotations, Fourier inte- 
grals are used to superimpose the various sinusoidal forces into step forces; that is. 



(A16) 

^oo 

= \ q.(t) dt 

(A17) 
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where q^(t) is the generalized ith coordinate. For a step function of amplitude 

q. equation (A17) jrields a frequency spectrum of 
i,s’ 

F(w) = Si (A18) 

io) 

with amplitudes, as given by equation (A16), that is, 
dqj = ^ F(w) e^"* do) 


or 

qj „ ^iwt 

dq. = — du) (A19) 

^ 2 ti ici) 

Equation (A12) yields the aerodynamic forces for any frequency. Hence, substitut- 
ing equation (A19) into equation (A12) and integrating over the whole spectrum, which 
constitutes a step variation in amplitude yields 



(A20) 


where the subscript s refers to step variation. Equation (A20) requires the evaluation 
of integrals of the type 


pOO 


0,2 e‘“* 

2-niu) 


dct) 


2irio) J 


(A21) 


The fimctions Lj and Mj 
ing quadratic expression in 


can be represented in their most general form by the follow- 
1/k ^illustrated for Ljj: 


Lj = f^O "*■ dQ + aj + (F + iG) dj — i- + (F i- iG) ^ (A22) 

L JL 


where the a-'s and d-j's are constants and (F + iG) represents the Theodorsen function 
J J 

which varies with k. Hence, the evaluation of equations (A21) requires the evaluation of 
the following basic integrals: 
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3 


') “ I. 


27rio)(k) 

(F + iG)fa;2 du) 
2irico(k)^ 




(A23) 


Remembering that 


a>=Xk 


and letting 


s = 


Vt 


leads to 


o)t = ks 


Thus, 


yj2 p«> ^2 iks 
Ij = Qr s^ ^ ^ 

J ^ 9 ^ J 1^2 J.oO 


27Tik(k)^ 


(j = 0, 1, 2) 


(A24) 


3 ’ ^ >-2 J_=o 2irik(k)^ 


(3 = 0 , 1 , 2 ) 


(A25) 


By remembering also that 


I 


«> -iks 

I— dk = H(s) 

_oo 2irik 


where H(s) is the imit step function, and that 


C (F + iG) _ e^^ ^ = 0(s) 
J.oo 2irik 
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where 0(s) is the Wagner function, equations (A24) and (A25) yield the following 
relationships: 


Iq = -Mr,s 

Jq = -do Hs) 

V • ^ 

Jl = -idi - cp{s) Qj. g 

T - 

^2 = ^ 2 ^ qr,s 


(A26) 


where q^, ^ and ^ represent, respectively, step-type variations in velocity and 
acceleration of the rth generalized coordinate. 


In the following equations, Ljq^. and 


Mjqj. 


are presented (as taken from refs. 22 


and 23) together with the respective results of step integration, denoted by ^oj^Ljqpj 


and 


pMjq,) 


Lfqr - 
[w^Liq. 


1 - i ^(F + iG) 

JK 


= -q 


r,s 


2 0(s) ^ q 


b ^r,s 


(A27) 


L2Qr = <1 - |[1 + 2(F + iG)] - ^(F + iG)> q^ 


(w2L2q 


= id 

r/g 2 


^(l + 2cp(s)) 

D 


\s - ^ %,s 


(A28) 


Ls^r = 


_i + i_i - i_ii 

7T k 77 k 77 


(F + iG) - AZlO(F + iG) 


2 7T 




r2 T 


T’’ b \ 77 77 




— <P(s) q 

77 ^ ^r,s 


(A29) 
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(A37) 
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s 4ir ^r,s 


(A38) 


Substituting equations (A27) to (A38) into the second term of equation (A20) yields 
the contribution of the controls to the generalized step forces. This relationship can be 
written in the form 


= upb^y^Uis) [Ai] + [AgjV + TTpb^vLs) [Bi] + [BgjV % + irpb^[E]i 


where 


-2 + 2 ^ 

IT > 




2t 

IT ^106 


^1 '^106 
b IT 


M = 


”^4 + '^10 


”^46 + 


[Bl] = 


^ ^ h ~ ^4 

b 7T 7T b 


o . ^2 " ^4 

^ -I 4- 

b \ b 7T 7T b ; 


'^116 _ ^^16 ^3 

7T 7Tb 


^i Fiia _ ^*^16 h 

b \ 7T 7T b , 


B2] = 


B2(1,1) B2(1,2) 

B2(2,1) B2(2,2) 
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where 


B 2 ( 1 , 2 ) 

B2(2,1) 

B2(2,2) 





+ 1 + 


2 p + T 4 
IT 


^5 h ' ^4 

7T b 


^1 "^46 ^ '^46 ^ *^56 h 

b 7T 7T 7T b 


[E] 


E(l,l) E(l,2) 
E(2,l) E(2,2) 


where 






E(l,2) = 

' ’ ' 7T TT b 


E(2,l) =4 
b 


_ 1 
b 2 


'^1 ^ ^2 
7T TT b 


- a, 


2 b 8 7T 


2/ TT 


^2 

4tt 


(A44) 




TT \ 0 2 / TT 4ir b 
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<1 


Figure 1.- General view and dimensions of Arava STOL transport. 






Arava fuselage center line 


Boom center line 



1 10 98 7654 321 


I I I ^ 

II 

Westwind fuselage center line 




Figure 3.- Strip allocations along wing and horizontal tail of 
Arava and Westwind aircraft. 









Aerodynamic forces 
Fj^^ = F^+F| Inertia forces 

Control forces 



Figure 6.- Distribution of various forces acting along' wing and locations of their 
resultants illustrating promotion of bending moment due to rigid-body— 
control-force interaction. 
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Figure 7.- Promotion of bending moment due to changes in inertia forces resulting 
from control forces acting at root of wing. 
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c 

Figure 8.- Variation with control rotation parameter of bending- moment 

Vg^Vc 

alleviation ratio at station 5 and maximum control rotation ratio. Arava 
transport with a single L.E.-T.E. active control system located at various 
strips along wing. 
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Without controls 

With controls 



Time, s Time, s 


Figure 11.- Variation with time of linear acceleration at center of gravity due to a (1 - cos) upgust. Westwind trans- 
port with a single L.E.-T.E. active control system located at various strips along wing and with = 0.5 rad. 









Figure 12.- Variation with control rotation parameter ^ — of acceleration alleviation 

VgjVc 

ratio at center of gravity and maximum control rotation ratio. Arava transport with a 
single L.E.-T.E. active control system located at various strips along wing. 
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Uc 

Figure 13.- Variation with control rotation parameter — of acceleration alleviation 

Vg/Vc 

ratio at center of gravity and maximum control rotation ratio. Westwind transport 
with a single L.E.-T.E. active control system located at various strips along wing. 
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Figure 15.- Variation with time of linear acceleration at point 1 and point 2 due to a 
(1 - cos) upgust. Westwind transport with a single L.E.-T.E. active control 
system located at various strips along wing and with Smax = 0-5 rad. 
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Figure 15.- Concluded, 
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Figure 16.- Variation with time of rotational acceleration at center of gravity due to a (1 - cos) upgust. Arava trans- 
port with a single L.E.-T.E. active control system located at various strips along wing and with 6max = 0*5 rad. 
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Figure 17.- Variation with time of rotational acceleration at center of gravity due to a (1 - cos) upgust. 
Westwind transport with a single L.E.-T.E. active control system located at various strips along 
wing and with 6max = 0-5 rad. 
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Figure 18 . r Variation with control rotation parameter — of flutter speed increase 

Vg/Vc 

ratio and maximum control rotation ratio. Arava transport with a single L.E.-T.E". 
active control system located at various strips along wing. 
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Figure 19.- Variation with control rotation parameter ^ — of flutter speed increase 

Vg/Ve 

ratio and maximum control rotation ratio. Westwind transport with a single L.E.-T.E. 
active control system located at various strips along wing. 
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Figure 21.- Variations with time of Awing bending moment at station 5 due to a (1 - cos) 
upgust. Arava transport with a single L.E.-T.E. active control system (using 
extended control law) located at various strips along Aving and Avith 6max = 0.5 rad. 
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Figure 22.- Variation with time of wing bending moment at station 10 due to a (1 - cos) 
upgust. Westwind transport with a single L.E.-T.E. active control system (using 
extended control law) located at various strips along wing and with S^ax = 0.5 rad. 
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Figure 23.- Variation with control rotation parameter ^ — of bending- moment 

Vg/Vc 

alleviation ratio at station 5 and maximum control rotation ratio. Arava trans- 
port with a single L.E.-T.E. active control system (using extended control law) 
located at various strips along wing. 
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Figure 24.- Variation with control rotation parameter 


Va/Vf 


of bending- moment 


alleviation ratio at station 10 and maximum control rotation ratio. Westwind 
transport with a single L.E.-T.E. active system (using extended control law) 
located at various strips along wing. 
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Figure 25.- Variation with time of linear acceleration at center of gravity due to a 
(1 - cos) upgust. Arava transport with a single L.E.-T.E. active control system 
(using extended control law) located at various strips along wing and with 
®max = 0*5 rad. 
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Figure 26.- Variation with time of linear acceleration at center of gravity due to a 
(1 - cos) upgust. Westwind transport with a single L.E.-T.E. active control 
system (using extended control law) located at various strips along wing and 
with 6 T-rk o-v = 0.5 rad. 






6c 

Figure 27.- Variation with control rotation parameter — of acceleration alleviation 

Vg/Vc 

ratio at center of gravity and maximum control rotation ratio. Arava transport with a 
single L.E.-T.E. active control system (using extended control law) located at 
various strips along wing. 
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ratio at center of gravity and maximum control rotation ratio. Westwind transport 
with a single L.E.-T.E. active control system (using extended control law) located 
at various strips along wing. 
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Figure 30.- Variation with time of linear acceleration at point 1 and point 2 due to a 
(1 - cos) upgust. Westwind transport with a single L.E.-T.E. active control 
system (using extended control law) located at various strips along wing and with 
^max ~ 
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Figure 31.- Variation with control rotation parameter 


of flutter speed increase 


ratio and maximum control rotation ratio. Arava transport with a single L.E.-T.E. 
active control system (using extended control law) located at various strips along 


wing. 
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Figure 32.- Variation with control rotation parameter — of flutter speed increase 

Vg^Vc 

ratio and maximum control rotation ratio. Westwind transport with a single L.E.-T.E. 
active control system (using extended control law) located at various strips along wing. 


88 



without controls 

With controls 






Figure 33.- Variation with time of various gust response parameters due to a (1 - cos) 
upgust. Arava transport with a single L.E.-T.E. active control system (using 
extended control law) located at strip 11, near center line of horizontal tail and 
with 6niax = 0-25 rad. 
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Figure 34.- Variation with time of various gust response parameters due to a (1 - cos) 
upgust. Westwind transport with a single L.E.-T.E. active control system (using an 
extended control law) located at strip 11 (at root of horizontal tail) and with 
5max = 0-5rad. 
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Figure 35.- Variation with control rotation parameter f — of various gust -alleviation 

Vg/Ve 

parameter ratios and maximvim control rotation ratio due to a (1 - cos) upgust. Arava 
transport with a single L.E.-T.E. active control system (using an extended control law) 
located at strip 11 near center line of horizontal tail. 




91 






Figure 36.- Variation with control rotation parameter j — of various gust -alleviation 

Vg/Vc 

parameter ratios and maximum control rotation ratio due to a (1 - cos) upgust. West- 


wind transport with a single L.E.-T.E. active control system (using an extended con- 


trol law) located at strip 11 at root of horizontal tail. 
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